Transport across a normal-superconducting interface: a novel probe of 
electron-electron interactions in the normal metal. 
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In this letter, we examine the effect of Coulomb interactions in the normal region of a normal- 
superconducting (N/S) mesoscopic structure, here the change from an attractive to a repulsive 
coulombic interaction, at the N/S interface, causes a shift in the order parameter phase. We show 
that this shift has a pronounced effect on Andreev bound states and demonstrate that the effect on 
Andreev scattering of non-zero order-parameter tails, can be used to probe the sign of the interaction 
in the normal region. 



Recent advances in material technology have enabled 
the fabrication of normal/superconducting (N/S) meso- 
scopic hybrid structures with well defined dimensions 
and interfaces [Q-Q. Due to the proximity of the nor- 
mal material to the superconductor the pairing field 
f(x) = {t/ji{x)ipi{x)) , in the normal region, decays to 
zero on the scale of a coherence length £ || . During the 
past decade this proximity effect has been extensively in- 
vestigated both experimentally and theoretically (see for 
example 

In contrast to the pairing field f(x), the effective 
electron-electron interaction, V(x), changes abruptly at 
the S/N interface, from an attractive interaction in the 
superconductor to either zero, a much diminished attrac- 
tive interaction or to a repulsive interaction, in the nor- 
mal (N) material. Consequently the order parameter, 
A = V(x)f(x), of a s-wave superconductor also changes 
abruptly at the S/N interface, as shown in figure [lj To 
date theoretical research into the transport properties of 
N/S interfaces have mainly considered the order param- 
eter in the normal region to be zero (figure |l]a) jl0],[llj . 

In this letter we consider the effect of an attractive 
or repulsive electron-electron interaction (V(x) ^ 0), as 
shown in [j]b and [j]c. The difference between figures [l]b 
and 0c is the tt phase shift in A(x), induced by the 
cross-over from an attractive to a repulsive interaction 
at the N/S interface. In what follows we examine how 
this phase shift affects the transport properties associ- 
ated with Andreev bound states. 

To investigate this regime we adopt a general scatter- 
ing approach to dc transport, which was initially devel- 
oped to describe phase-coherent transport in dirty meso- 
scopic superconductors |12|. For simplicity in this let- 



ter, we focus solely on the zero- voltage, zero-temperature 
conductance, for the structure shown in figure |^. In 
the linear-response limit, at zero temperature, the con- 
ductance of a phase-coherent structure may be calcu- 
lated from the fundamental current voltage relationship 
©0, 



mal reservoir i to the voltage differences (vj — v), where 
v = ji/e and the sum is over the 2 normal leads con- 
nected to the scattering region. The 's are linear com- 
binations of the normal and Andreev scattering coeffi- 
cients and in the absence of superconductivity satisfy 
2i=i a ij = Sj=i a ij — m which case the left hand 
side of equation pi becomes independent of v. In units of 
2e 2 /h , ati = Ni + Rf - R° and a mi =T$-T§, 

where T,4, TS are Andreev and normal transmission co- 
efficients from probe j to probe i, Rf, i?P are Andreev 
and normal reflection coefficients from probe i and iVj is 
the number of open scattering channels in lead i. 

Setting I\ = —1-2 = I and solving equation |l| for the 2 
probe conductance yields 



G = 



(Vi - V 2 ) 



= T 2 \ 



2{RfRf - T£T{ 



Rf 



Rf 



(2) 



where G is the conductance in units of -=f-. As noted 
in pa ] the various transmission and reflection coeffcients 
can be computed by solving the Bogoliubov - de Gcnncs 
equation on a tight-binding lattice of sites, where each 
site is labelled by an index i and possess a particle (hole) 
degree of freedom tp(i) (</>(«)). In the presence of local 
s-wave pairing described by a superconducting order pa- 
rameter Aj, this takes the form, 



Eipi = Eiipi - J2s 7 (V-i+a + ipi-s) + A^i 
Efa = -e 4 4 + (<t>i+6 + <Pi-s) + A*ip u 



(3) 



In what follows, the Hamiltonian of eq.(g) is used to 
describe the structure of figure ^, where for all i, the on- 
site energy ej = cq. In the S-region, the order parameter 
A, is set to a constant, A, = Ao, while in the normal 
region Aj is approximated by, 



Aj = (tanh(j - L n ) + 1) . 



(4) 



j(vj - v), 



(1) 



The above expression relates the current Ii from a nor- 



The nearest neighbour hopping element 7 merely fixes 
the energy scale (ie the band- width), whereas eo deter- 
mines the band-filling, and L n is the length of the normal 
region. In what follows we choose 7 = 1 and A = 0.1. 
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By numerically solving for the scattering matrix of equa- 
tion ||, exact results for the dc conductance can be ob- 
tained and therefore the effects of a repulsive/attractive 
coulomb interaction in the normal region can be exam- 
ined. 

For simplicity in this letter we consider the transport 
properties of the structure shown in figure |^. Consisting 
of two normal, semi-infinite, crystalline leads, 20 sites 
wide, joined by a scattering region. The normal scat- 
tering region is 40 sites long and the superconducting 
region is of length L$- To form Andreev bound states 
we create quasiparticle confinement in the area in front 
of the superconductor by introducing a weak point-like 
contact between the left lead and the scattering region, 
and ensuring that there is no or very little quasiparti- 
cle transmission through the superconductor. For this 
reason the superconductor length was set to, L$ = 150 
sites. However care must be taken since the decay length 
of sub-gap states in the superconductor increases with 
quasiparticle energy. At energies close to the gap energy 
the decay length is long enough for transmission through 
the superconductor, see figure [| The weak contact is 
created by placing a potential barrier, one site wide, be- 
tween the left lead and the normal scattering region, via 
a mean potential U added to the on site energy €q. Fig- 
ure H shows a plot of the transmission coefficient of the 
barrier as a function of U . From this plot we see that 
a for low quasiparticle transmission through the barrier 
a high potential is needed. For these reasons the barrier 
potential was set at U = 20. 

These conditions produce discrete states within the 
normal regions. Allowing the formation of Andreev 
bound states, when these energy levels become populated 
reasonances appear in the conductance of the system, 
analogous to Breit-Wigner resonances ||l7| . By plotting 
G as a function of quasiparticle energy (all energies are 
with respect to the Fermi energy) we are able to inves- 
tigate transport resonances in the conductance due to 
the formation of Andreev bound states. Figure || shows 
the conductance function of quasiparticle energy 

for the three proximity tails corresponding to the van- 
ishing, attractive and repulsive electron-electron interac- 
tions respectively. These represent the central result of 
this letter. We see that the introduction of a proxim- 
ity tail has the effect of shifting the energy at which the 
Andreev bound states occur. A positive proximity tail 
causes a shift in the resonance energies to the right of 
the zero-tail spectrum shown in figure ra, whereas a nega- 



tive proximity tail shifts the spectrum to the left. These 
results suggest a novel method for detecting the sign of 
the electron-electron interaction in the N-metal. Namely 
by suppressing the order parameter in the normal region, 
the resonances will either shift to higher energies indi- 
cating a repulsive interaction, or lower energies in the 
case of an attractive interaction. This suppression can 
be achieved by, for example applying a magnetic field or 
via a control current in the superconductor. 

This work is funded by the EPSRC and the E.U. TMR 
programme. 
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FIG. 1. a) Self-consistent order parameter profile, V(x)=0 
in the normal region, b) Self-consistent order parameter 
profile, V(x) and f(x) both finite in normal region, c) 
Self-consistent order parameter profile including Coulomb in- 
teractions. 
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FIG. 2. A N-S system with a tunnel barrier present. Such 
a system will have Andreev bound states present in the region 
between the interface and the tunnel barrier. 
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FIG. 3. Graph showing the transmission coeffcient as a 
function of Superconductor length for various quasiparticle 
energies. No Barrier or proximity tail are present. Structure 
width is 20 sites. In this graph it can be seen that at high 
quasiparticle energies transmission will occur through the su- 
perconductor. 



FIG. 4. Graph showing the effect of barrier potential. The 
barrier is one site wide, structure width is 20 sites. The graph 
shows transmission as a function of barrier potential for a 
quasiparticle energy, E=0.00. 
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Effect Df j n I jtatio] in the Dtde] Parai jtei Phase n tl 3 Nori lal Metal 
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FIG. 5. Graphs of Conductance as a function of quasipar- 
ticle energy for different phases of the order parameter tail in 
the Normal metal. 
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